Accepted for Publication in ApJ Letters 

Preprint typeset using LAT^X style cmulateapj v. 5/2/11 



THE LACK OF TORUS EMISSION FROM BL LACERTAE OBJECTS: AN INFRARED VIEW OF 

UNIFICATION WITH WISE 
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ABSTRACT 

We use data from the Wide-Field Infrared Survey Explorer ( WISE) to perform a statistical study 
on the mid-infrared (IR) properties of a large number (^10 2 ) of BL Lac objects — low-luminosity 
Active Galactic Nuclei (AGN) with a jet beamed toward the Earth. As expected, many BL Lac 
objects are so highly beamed that their jet synchrotron emission dominates their IR spectral energy 
distributions. In other BL Lac objects, however, the jet is not strong enough to completely dilute 
the rest of the AGN emission. We do not see observational signatures of the dusty torus from these 
weakly beamed BL Lac objects. The lack of observable torus emission is consistent with suggestions 
that BL Lac objects are fed by radiatively inefficient accretion disks. Implications for the "nature 
vs. nurture" debate for FR I and FR II radio galaxies are briefly discussed. Our study supports the 
notion that, beyond orientation, accretion rate plays an important role in AGN unification. 
Subject headings: accretion, accretion disks — BL Lacertae objects: general — infrared: galaxies 



1. INTRODUCTION 

The standard orientation-based unification paradigm 
posits that every Active Galactic Nucleus (AGN) is com- 
prised of the same basic components: an accretion disk 
around a supermassive black hole, a broad emission line 
region (BELR), an obscuring torus, and a narrow line 
region; radio- loud AGN also launch large-scale relativis - 
tic jets (e.g., lAntonuccl 119931 : lUrrv fe Padovanl 119951 ). 
The torus is a key component in this unified model, as 
it can block the BELR along certain lines of sight, and 
its dust also reprocesses UV/X-ray radiation into the in- 
frared (IR). In simple unification, the only other way to 
hide the BELR is also via a geometric argument: emis- 
sion lines are outshown by a relativistic jet beamed to- 
ward Earth. 

In addition to orientation, intrinsic differences among 
AGN also play a r ole in AGN unification. For example, 
lUrrv fc Padova ni (1995) conclude their seminal review 
with 10 outstanding questi ons, including one rega rding 
Fanaroff- Riley galaxies (FR; Fanaroff fc Rile ylfl 9741 ): "do 
FR Is have broad emission line regions?" The significance 
of this question is that if FR I galaxies lack BELRs, 
then intrinsic properties might play a prominent role in 
driving the so-called FR I /II dichotomy (in which the 
more powerful FR II galaxies have edge-brightened radio 
lobes). Indications so far are that emission line lumi- 
nosities from FR I galaxies are 5-30 times weaker than 
for FR II galaxies, and FR I galaxies may also have 
weaker tori (|Zirbel fc Baumlll99alChiaberge et al.lll999t 
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Donato et al.l 12004 ; alth ough also see iCao fc Rawlingsl 
2004HLeipski et alll200ll . However, the requisite obser- 



vati ons and their interpretation are techn ically challeng- 
ing (jCapetti et al.ll2000t lE^ans et al.ll2006l ). and there are 
still many ongoing "nature vs. nurture" debates. 

We can turn to beamed versions of FR I and FR II 
galaxies for additional guidance — i.e., BL Lac objects 
and flat spectrum radio quasars (FSRQs), respectively, 
which are collectively called blazars. Broad emission 
lines and sometimes dusty tori are routinely detected 
in the spectral energy distributions (SEDs) of FS RQs, 
but less often from BL L ac objects (e.g. lAbdo et all 
l2010t IGiommi et al.l 1201 lh . However, it is generally 
still not observationally clear if BL Lac BELRs and 
tori appea r weak solely due to dilution by je t emis - 
sion (e.g., IChen fc SharJ WHi IMalmrose etail [20T1 . 
or rather if they r e ally are intrinsically weak (e.g., 
iGhisellini et alJl20lH fSbarrato et al.ll201l( i. Constraints 
on blazar BELRs and tori are furthermore important, as 
both components are potential sources of seed photons 
for Comptonized gamma-ray emission, especially for FS- 
RQs and low synchrotron peaked (LSP) BL Lac objects. 
Note, gamma-ray emission from high synchrotron peaked 
(HSP) BL Lac objects is g enerally consistent with syn- 
chrotron self Compton (e.g. lAbdo et al.ll2010D . 

The preliminary data rel ease of the Wide-fiel d Infrared 
Survey Explorer (WISE; Wrig ht et al.l 120 lOl) recently 
opened a new multiwavelength window over a large area 
of the sky. In this letter, we investigate the tori of low- 
luminosity radio galaxies by examining BL Lac objects 
detected by WISE. We describe our BL Lac sample in 
Sj2] In [J3] we describe their IR properties, and results 
are discussed in £|4j All spectral indices are defined as 
f v oc v~ av , and we define the broad-band radio to IR 
spectral index a r i — — log(L r /Lj)/4.25, where L r , and 
Li are monochromatic luminosities at 5 GHz and 3.4 /im 
rest-frame, respectively. We adopt the following cosmol- 
ogy: H = 71 km s" 1 Mpc" 1 ; Cl M = 0.27; ft A = 0.73. 
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Fig. 1. — Sample SDSS optical spectra for four BL Lac objects (left), and corresponding SEDs (right). For illustrative purposes, we 
assume z = 0,2 for the top two rows. WISE data points are shown as open squares, and the other SEP multiwav clcngth data are assembled 
as in §4.3.5 of I Plotkin etaLl H2011al) . For comparison, we include a typical quasar SED (from Richards ct al. 2006, R06), which shows 
thermal torus dust emission peaking redward of 1—2 (im (3.0— 1.5xl0 14 Hz). 



2. THE BL LAC SAMPLE 

We start with 59 confirmed radio-lou d BL Lac objec ts 
from the catalog of iPlotkin et~all (|2010l . hereafter [Pl3) Q 
which were selec ted from Sloan D igital Sky Survey spec- 
troscopy (SPSS: rYo7ket~a^l2000l ). iPlOl require no SDSS 
spectrum to show an emission feature with REW > 5 A, 
and the Ca II H/K break, C, must be small er than 40% 
(e.g. IStocke et al.l[l99ll IMarcha et aflfl996l) . The H/K 
break (i.e., the fractional change of continuum flux sur- 
rounding 4000 A rest-frame) quantifies the host galaxy 
(HG) contribution to each SDSS spectrum. A spectrum 
with a more highly beamed jet will show smaller C be- 
cause th e jet is brighter com pared to the (unbeamed) HG 
(see, e.g. lLandt etal|[2002l ). Of the ~60% of BL Lac ob- 
jects for which [P 101 can determine redshifts, most (and 
essentially all at z < 0.5) are derived via absorption lines 
from H G sta rlight and are thus among the most weakly 
beamed |P10l BL Lac objects (see Figure [T]). 

6 A subset of radio-quiet B L Lac candidat es from IP101 with 
WISE coverage are discussed in lWu et afl pOllT l. IWu et alj pmT l 
show that many are likely low-redshift analogs to weak line quasars 
(WLQs), which are high-redshift (z > 3) quasars wit h unusually 
weak or missing BELRs but normal dusty tori (see ILane et al.l 

120m 



We correlate the IP10I BL Lac objects to the prelimi- 
nary WISE data release using a 3" search radius. We re- 
quire WISE detections with S/N > 3 in the Wl[3.4 //m], 
W2[4.6 //m], and W3[12 //m] bands. Requiring S/N > 3 
for the W4[22 fim] band would reduce our sample by at 
least a factor of two. We retain 157 objects (see Fig- 
ure [TJ . To create comparison samples of normal quasars 
and early-type galaxies (ETGs; BL Lac objects live per- 
haps universally in large ellipticals. lUrrv et al.ll2000h . we 
also correlate 10 5,783 spectroscopically confirmed Type I 
SDSS quasars ()Schneider et al.1 12010T) and 8 666 (inac- 
tive) early SPSS galaxies (jBernardi et al.ll2003[ ) to WISE 
applying the same S/N constraint as above. All IP10I 
BL Lac objects have z < 2, so we also restrict the com- 
parison quasars to z < 2. We further consider only 
quasars with central black hole masses 10 8 < M^/Mq < 
10 9 (masses from lShen etaI]|20llT ), and ETGs with stel- 
lar velocity dispersions 200 < (Jdisp < 320 km s -1 . The 
abo ve values of Mbh and adisv are typical for BL Lac HGs 
(see lLeon-Tavares et al.l mm IPlotkin et al.ll2011b[ ). We 
are left with 13,881 quasars and 747 ETGs. Throughout, 
we convert WISE magnitudes to flux densities (/„) using 
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Fig. 2.— Distributions of W1-W2 (a) and W2-W3 (b) for lPlOl 
objects lacking redshift s. Pa nels (c) and (d) show redshift vs. IR 
colors for the remaining P1Q BL Lac objects (blue squares). Open 
symbols indicate less reliable redshifts. We also show the IR colors 
of our comparison quasars (contours) and ETGs (filled circles). 

the the published (iso) zero pointsQ 

3. BL LAC INFRARED PROPERTIES 



The WISE colors, W1-W2 and W2-W3, for the IPlOl 
sample and for our comparison quasars and elliptical 
galaxies are shown in Figure [5] as a function of redshift. 
The dusty torus of quasars makes their IR colors red- 
der than for ETGs. The first-order dependence of the 
quasar IR colors on redshift is simply due to the torus be- 
ing redshifted through the WISE filters. Although some 
BL Lac objects have WISE colors similar to quasars, 
BL Lac objects tend to populate the bluer edge of the 
quasars' color space. For most lower- redshift BL Lac 
objects, this is because of HG contamination (elliptical 
galaxy SEDs emit most of their radiation in the near-IR 
and appear blue in the mid-IR; see Figure [TJi). However, 
beamed synchrotron emission should completely domi- 
nate the IR emission from the hi gher redshift objects 
and those lacking re dshifts (see, e.g. lPadovani et aD2006t 
iChen fc Shanl I2011I ) . We thus expect any overlap be- 
tween the IR colors of highly beamed BL Lac objects and 
SDSS quasars to be coincidental, since different mecha- 
nisms (i.e., synchrotron vs. thermal dust emission) are 
producing their IR flux. 

3.1. Constraints on BL Lac Dusty Tori 

BL Lac objects bridge the IR color space between nor- 
mal galaxies and quasars (Figure^)- One might expect 

7 http: / / wise2.ipac.caltech.edu / docs / release/prelim/expsup / figures / 
sec4_3gt4.gif. 



the most weakly beamed BL Lac objects to show contri- 
butions from both HG starlight and from thermal dust 
emission from the torus. Then, the most weakly beamed 
BL Lac objects would look more similar to quasars that 
appear extended in their SDSS images (red plus signs). 
However, the weakly beamed BL Lac objects' IR colors 
are too blue for a significant contribution from the torus. 
In this section we determine if the lack of torus signa- 
tures is simply because even weakly beamed jets (and 
HGs) can still outshine the torus. 

All ETGs in our comparison sample have z < 0.3, so 
we restrict ourselves to 28 BL Lac objects with z < 0.3 
in this section. All 28 objects appear as point sources in 
their WISE images. There is thus little concern that the 
WISE photometry includes different fractions of the total 
HG flux as a function of redshift. We also verify that 
they show no significant correlation between redshift and 
either WISE color. All 28 objects have z > 0.1, so, for 
uniform comparison, we only consider ETGs and quasars 
from 0.1 < z < 0.3 (530 and 470 objects, respectively). 

There are three emission components that must be 
modeled in the IR: 1) starlight from the HG; 2) beamed 
jet synchrotron emission; and 3) dusty torus emission. 
To test if the dusty torus is present, we run Monte Carlo 
(MC) simulations. These simulations involve building 
a set of probability density functions (pdfs; see H3.1.1I 

t j3.1.3[) . from which we take 10 6 random draws to 
estimate each IR-emitting component's monochromatic 
luminosity in each WISE filter. From these luminosi- 
ties we synthesize the expected WISE color space occu- 
pied by BL Lac populations with and without the torus. 
Throughout we estimate all luminosities at z = 0.2 
to minimize K-correction uncertainties, and we restrict 
comparison samples to a narrow redshift range near 
z = 0.2. For all pdfs modeled as (log)-normal distribu- 
tions below, we first check that their observed parameter 
distributions are indeed approximately Gaussian. 

3.1.1. IR Emission from Host Galaxy Starlight 

We consider the observed properties of 164 ETGs with 
z = 0.22 ± 0.05 (the median of our low-redshift BL Lac 
sample). We first build a pdf for the monochromatic 
luminosity in the W2 filter, assuming a log-normal dis- 
tribution with {logL v ) hg w2 = 29.82±0.14ergs" 1 Hz" 1 . 
Then we make pdfs for HG colors using normal dis- 
tributions with (Wl - W2) = 0.30 ± 0.066 mag and 
(W2 - W3) = 1.58 ± 0.43 mag. Randomly drawing from 
each distribution allows us to estimate HG monochro- 
matic luminosities in each WISE filter. This method 
implicitly accounts for measurement uncertainties. 

3.1.2. Jet Synchrotron in the IR 

We begin with an unbeamed 5 GHz radio core luminos- 
ity (vL v )- r , and we use the broad-band radio-IR spectral 
index a r i to estimate an unbeamed IR jet luminosity in 
the Wl filter, {yL v )^ wl . Beamed IR jet luminosity is 

then calculated as [yL v ^- wl = (yL v )- wl S 2+awis ^ (see 

ILind fc BlandfOTdl[T985l) . where 5 = [7 (1 - ficose)}' 1 is 
the Doppler parameter, 7 = (1 — /3 2 )~ - 5 , 6 is the viewing 
angle (8 = corresponds to a perfectly aligned jet), and 
3 is the jet velocity normalized to the speed of light. We 
then use (vL^)^ wl and the local IR spectral index a W i se 
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Fig. 3. — (a) Color-color diagram for 28 low-rcdshift BL Lac objects (blue squares), 150 comparison quasars with SDSS point-like 
morphologies (circles), 320 quasars that appear extended in SDSS imaging, and 530 ETGs (crosses). In the remaining panels, contours 
show expected BL Lac colors for MC simulations including only a beamed jet (b), a jet and host galaxy (c), and a jet, host galaxy, and 
torus (d) (see text for details). Quasars and ETGs are omitted from panels (b)-(d) for clarity. 



to estimate beamed jet luminosities in the W2 and W3 
nlters0 

To estimate beamed jet luminosities for a random view- 
ing angle, Lorentz factor, and jet power, we randomly 
draw from the following pdfs: a log-normal distribution 
for {vL v ) jr with {\ogvL v ) jr = 40.77 ± 0.69 ergs" 1 (de - 
rived from Equation 2 and Table 4 of IWu et all 120071 ). 
a normal distribution f or 7 with (7) = 7 ± 0.7 (e.g., see 
iMerloni fc Heind 120071 ) , and a uniformly distributed 
from 0-40°. For a wise , we fit power laws to the Wl, W2 
and W3 flux densities (in log f u — log v) for the 69 IPIOI 
BL Lac objects lacking redshifts. Based on those mea- 
surements, we then randomly draw a W i se from a nor- 
mal distribution with (a W i se ) = 0.60 ± 0.43. We simi- 
larly measure a r i for those 69 BL Lac o bjects (assuming 
z = 0.2) using the radio fluxes in IPIOI and the Wl fil- 
ter. We then randomly draw a r i from a normal distribu- 
tion with (a r i) = 0.40 ± 0.11. Finally, once we estimate 
luminosities in each filter as described in the previous 
paragraph, we add random noise to each luminosity by 
assuming a L /L = ±0.05, 0.05, and 0.1 in the Wl, W2, 
and W3 filters, respectively (based on typical WISE flux 
measurement uncertainties for our full BL Lac sample). 

8 We assume a w i s< , and a r ; are not strongly affected by beaming, 
which is not strictly true but a reasonable approximation once the 
jet is highly enough beamed to appear as a BL Lac object. 



3.1.3. IR Emission from a "Typical" Obscuring Torus 

We employ a realistic coupling between the jet and ac- 
cretion flow to estimate IR luminosities from the dusty 
torus. First we estimate X-ray luminosity L x (from 2— 
10 keV) usi ng (vL^)- r from a bove and the following re- 
lations from IMerloni Sz Heinz! ([20071 ). based on a sample 
of low-luminosity AGN with mechanical jet power esti- 
mates: 



log 



Li 



^Edd 



Hog 




_^Edd_ 





log [5/0. 2; 



(0.78 ±0.36) 



and 



logL fem = (0.81 ±0.11) log (uL u ). +11.9 



+4.1 



(1) 



(2) 



where L^m is the kinetic power of the jet in erg s , 
L Edd = 1.3 x 1O 38 (M/M ) erg s" 1 is the Eddington 
luminosity, and B is fraction of the bolometric lumi- 
nosity emitted from 2-10 keV. We note that Equa- 
tions [T] and [2] are consistent with the theoretical scalings 

for jet dominated accretion flows, HzaL,,^. r / L, 
L kin/L E dd) 17/12 and (L k in /L Edd ) oc {L x /L Edd )°- 5 (e.g. 



* (hi 



OC 



Falcke &: Biermarml I1995D . The above scalings are ap- 
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plicable to hard state X-ray binaries and their super- 
massive analogs, in cluding BL Lac objects (see, e.g., 
IPlotkin et al"1l2011alL 

We then assume X-ray luminosity is a good tracer for 
the t orus' IR luminosity, using the IR/X-ray correlation 
from lGandhi et all (pOOl : 



log (iW 

\ ±u / 12.3/im 



(1.11 ±0.07) log (i^) 
f (0.19 ±0.05), 



(3) 



where (^£„)i2.3^m is the IR luminosity at 12.3 /zm, and 
the X-ray luminosity, L x , is from 2-10 keV. 

Combining Equations[TJ[2j and[3]leaves us the following 
relation for the IR torus emission at 12.3 /im, assuming 
an unbeamed radio core luminosity at 5 GHz: 



V ±u / 12.3/jm 



1.84 log 



(vL„\ 

Vio 40 /. 



1.11 log 



B 
OA 



1.16 log 



0.43, 



^Edd 



10~ 6 from KS tests comparing their W1-W2 and W2- 
W3 distributions to the synthesized jet). 

To explain the SDSS population of weakly beamed BL 
Lac objects (blue squares), we favor the synthesized col- 
ors including the jet and HG (Figure |3t) over the colors 
that also include the torus (Figure |3jl) • The simulations 
including the torus predict that BL Lac colors should 
be weighted toward redder IR colors, while the simula- 
tions with only the jet and galaxy predict IR colors more 
consistent with observations. KS tests show that the 
W1-W2 and W2-W3 distributions for 28 BL Lac objects 
with 0.1 < z < 0.3 and the synthesized colors including 
the torus are statistically different (p ~ 10~ n and 10~ 8 , 
respectively), while similar parent distributions are not 
highly excluded between the observed colors and the syn- 
thesized colors without the torus (p ~ 0.05 and 0.02). 



4. DISCUSSION AND CONCLUSIONS 



where we assume B — 0.1 (to conservatively err on the 
side of underestimating the torus flux). Finally, we ex- 
trapolate torus luminosity to the W2 filter assuming an 
IR spectral index of 0.15 between 1 2.3 and 4.6 /xm (esti- 
mated at z = 0.2 from the average [Richards et al.l [20061 
quasar SED), and we estimate IR luminosities in the Wl 
and W3 filters using pdfs based on the WISE colors of 
150 quasars (0.1 < z < 0.3) with point-like SDSS mor- 
phologies. 

For the MC simulations, we build the following ad- 
ditional pdfs. For luEdd we assume BL Lac black hole 
masses are log-norm al with (logM) = 8.54 ± 0.40 M 
(|Plotkin et allboilblh Since Equation [3] relies on several 
assumptions that propagate non-linearly, we randomly 
add intrinsic scatter to each estimate of (log^L„) 12 3 
assuming a generous o~i n t = ±1.5 dex. For the pdfs 
for Wl — W2 and W2 — W3, we first randomly draw 
Wl — W2 from a normal distribution with (Wl — W2) — 
1.03 ± 0.10 mag. Since quasar WISE colors are cor- 
related, we then use the best-fit relation W2 — W3 = 
1.19 + 1A7(W1 - W2) to build the W2 - W3 pdf, and 
we add random noise using the root-mean-square (rms) 
scatter about the best-fit regression, a rms ± 0.23 mag. 

3.1.4. Lack of observational signatures from the dusty torus 

We perform 10 6 MC simulations to estimate luminosi- 
ties from each IR-emitting component in each WISE fil- 
ter. Figure [3] shows the synthesized WISE colors if only 
one, two, or all three components are included. Note, 
the simulated jets encompas s a range of viewing angles 
6 < 40° and Lorentz factors. iMassaro et al.l ([201 ID show 
that blazars fall in a distinct region of IR color space, 
which th ey ca ll the ' WISE blazar strip.' The most highly 
beamed IPIOI BL Lac objects (i.e., the 69 lacking red- 
shifts) also fall within the strip, and a KS test shows 
that the W1-W2 and W2-W3 distributions for those 69 
BL Lacs and our synthesized 'blazar strip' (panel b) are 
not significantly different (p ~ 0.4 and 0.6, respectively). 
However, for the 28 weakly beamed objects, a bluer IR 
component (i.e., the HG) is necessary (p ~ 10~ n and 



3.78 • 10 46 / We conclude from the above that our weakly beamed 
BL Lac objects lack observational signatures of the dusty 
torus in the mid-IR. Thus, BL Lac dusty tori appear 
weaker than the tori of luminous quasars (i.e., one can- 
1 not simply scale a normal quasar torus down to BL Lac 
luminosities). Another way to express this conclusion is 
that BL Lac tori have different properties than the tori 
of normal quasars. It is reasonable to extrapolate this re- 
sult to highly beamed BL Lac objects, if beaming is pri- 
marily a geometric argument. If one interprets the torus 
as an extension of the BELR (e.g., lElitzur fc Shlosmanl 
2006), then our study suggests that BL Lac BELRs are 
also intrinsically weak. We stress that this is a statisti- 
cal conclusion. The simulations including the torus do 
predict a small population of weakly beamed BL Lac ob- 
jects with WISE colors matching the observations, but 
that population is not large enough to explain the IR 
colors for the BL Lac population as a whole. However, 
we cannot exclude the presence of the torus/BELR from 
every single BL Lac object. Indeed, some BL Lac objects 
can show very weak broad emission lines in their optical 
spectra (especially some LSPs). 

One explanation for weaker BL Lac tori is that their 
physical properties (i.e., torus geometries, dust proper- 
ties, ionizing SEDs, etc.) are different than for luminous 
quasars. Another interpretation is the same one invoked 
for optically du ll AGN and "naked" Seyfert galaxies (e.g., 
lHawkinsl[200l iTrump et al.ll201lD . as follows. Consider 
a scenario where the BELR an d torus are fed by a radia - 
tively driven disk wind (e.g.. iMurray fc Chiang! [1998). 
At sufficiently low accretion rate (below a few percent 
L/LEdd), the inner region of the accretion disk is re- 
placed by a radiatively inefficient accretion flow (RIAF) . 
The RIAF is then no longer able to support a sufficient 
wind to populate the BELR and torus (e.g., iNicastrol 
[2000h . In this scenario, the strength of the BELR and 
torus is intimately connected to the efficiency of the ac- 
cretion flow and therefore accretion rate0 Indeed, such 
an "accretion rate dichotomy" has already been sug- 
gested to explain the BL La c/FSRQ divide (with BL 
Lac objects fed by RIA Fs, e.g.. lBottcher fc Dermerll2002l : 
iGhisellini et al.1 [2009( 1. Such a scheme may also be re- 

9 Note, WLQs and their weak BELRs require a different expla- 
nation b ecause WLQs sho w normal dusty tori and are radiatively 
efficient l|Lane et al.|[2uTJ) . 
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sponsible for the FR I/II divide (e.g.. lGhisellini fc Celottil 
120011: iWold et all 120071) . In fact, an important (sup- 
porting) point is that the classification of unbeamed ra- 
dio galaxies based on high-excitation vs. low-excitation 
emission lines (which are likely connected to the accre- 
tion flow) is more physically meaningful than the older 
FR I/II scheme based on radio morphology /jet powe r 
(see, e.g.. Uackson fc Wall|[l999t lHardcastle et all 12009). 
Furthermore, such a division is naturally expected from 
the hysteresis displayed by accreting stellar mass black 
holes, where changes in X-ray spectral states for individ- 
ual objects (i.e., at fixed orientations) may similarly be 
related to the efficiency of the inner accretion flow. 
WISE offers a new perspective on the strength of the 



obscuring torus in low-luminosity radio galaxies. The 
WISE view of BL Lac objects is consistent with other 
(independent) indications of an accretion rate dichotomy 
among AGN, adding a complementary piece to the puz- 
zle. Further constraints on such a division in accretion 
mode is an important step for better understanding AGN 
unification. 
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